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Abstract 

We study a recently introduced model [8,9] which consists of positive and negative parti- 
cles on a ring. The positive (negative) particles hop clockwise (counter-clockwise) with rate 
1 and oppositely charged particles may swap their positions with asymmetric rates q and 
1. In this paper we assume that a finite density of positively charged particles p and only 
one negative particle (which plays the role of an impurity) exist on the ring. It turns out 
that the canonical partition function of this model can be calculated exactly using Matrix 
Product Ansatz (MPA) formalism. In the limit of infinite system size and infinite number 
of positive particles, we can also derive exact expressions for the speed of the positive and 
negative particles which show a second order phase transition at q c = 2p. The density pro- 
file of the positive particles on the ring has a shock structure for q < q c and an exponential 
behaviour with correlation length £ for q > q c . It will be shown that the mean-field results 
become exact at q = 3 [8] and no phase transition occurs for q > 2. 

PACS number: 02.50.Ey, 05.70.Fh, 05.70.Ln 

Key words: Matrix Product Ansatz (MPA), Non-equilibrium phase transition, Shock profile 



e-mail: JAFAR@theory.ipm.ac.ir 



1 Introduction 



One-dimensional driven diffusive systems have provided a rich frame-work for the study of 
many interesting phenomena such as phase transitions, shock structures and the spontaneous 
breaking of translational invariance in the context of the non-equilibrium statistical mechanics 
[1-4]. Moreover, from the Mathematical point of view, some of them are related to the inte- 
grablc quantum chain Hamiltonians and can be mapped to the other non-equilibrium models, 
for instance, interface growth and traffic flows [5,6]. 

In this paper, we study a recently introduced exclusion model which exhibits both phase tran- 
sition and spatial condensation of particles [8,9]. In this model two types of particles, called 
positive and negative, occupy the sites of a periodic one-dimensional lattice of length L. The 
particles are subjected to hard-core exclusion, so that there are three possible states at each site: 
empty, occupied by a positive particle, or by a negative one. The positive (negative) particles 
hop to their immediate right (left) site with rate 1 provided that it is empty. Adjacent positive 
and negative particles also exchange their positions with asymmetric rates q and 1. Therefore, 
during the infinitesimal time step dt, any bond (i, i + 1) (with 1 < i < L) evolves as follows: 
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with 


rate 
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with 
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-> (-)(+) 


with 
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With periodic boundary conditions, the stationary state of this model has been studied by 
Monte-Carlo simulation, mean-field approximation, and recently using analytical approaches in 
the neutral case in which the density of the positive and negative particles are equal [8-11]. It 
has been shown that in the neutral case the model possesses two different phases depending on 
the reaction rate q [10], in contrast to the analytical results of [11] which indicate that there 
is no such a transition. This model has also been studied with open boundaries in two special 
limits [12]. 

Here we study this model in the charged case in which a finite density of positive particles p and 
only one negative particle are present on the ring. Using the Matrix Product Anstaz (MPA) 
formalism first introduced in [7], we can derive the exact analytic expression for the partition 
function of this model. Exact calculations show that in the limit of infinite system size and 
infinite number of positive particles the model has two different phases for q < 2, and the tran- 
sition point q c depends on the density of the positive particles. Exact expressions for the speed 
of the positive and negative particles can also be computed in this limit. We will see that these 
expressions have different behaviours depending on the value of q. For q < q c , the speed of the 
positive and negative particles, as a function of q, increase linearly with q from zero. However, 
for q > q c the speed of the positive particles remains constant while the speed of the negative 
particle still increases as a function of both p and q. The speed of particles as a function of 
q is continuous at the transition point although its derivative changes discontinuously at this 
point. Also for q > q c it turns out that the density profile of the positive particles has an 
exponential behaviour with a characteristic length £ = |ln^| _1 , while the system presents a 
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shock, i.e. a sharp discontinuity between a region of high density of particles and a region of 
low density for q < q c . At the transition point the correlation length diverges and the model 
shows a second order phase transition. We will show that at q = 3 the stationary probability 
of all possible configurations become equal. At this point no correlation exists and mean-field 
results are exact [8]. 

This paper is organized as follows. In the section 2, we obtain the exact expression for the 
canonical partition function of the model using the MPA. In section 3, we derive exact expres- 
sions for the speed of both kinds of particles, and also for the density profile of the positive 
particles on the ring. In the last section we compare our results with those obtained from the 
neutral case. 



2 Expression of the canonical partition function using the MPA 

According to the MPA formalism [7], the stationary probability distribution P({C}) of any 
configuration {C} can be expressed as a trace of a product of non-commuting operators. For 
the model proposed here we assume that there are M positive particles and only one negative 
particle on a ring on length L. Since this model is translationally invariant, we can always keep 
the single negative particle at site L and write the normalized stationary probability distribution 
as 

P({C}) = ^—Tr( ]J (nD + (1 - Tl )E)A) (2) 

where T{ = 1 if the site i is occupied by a positive particle and Tj = if it is empty. The 
matrices D, A and E which stand for the presence of a positive, negative particle, and a hole 
satisfy the following algebra introduced in [10] 

qDA-AD = D + A (3) 
DE = E (4) 
EA = E. (5) 

The normalization factor Zlmi which plays the role analogous to the partition function in 
equilibrium statistical mechanics, ensures that J2aii conf. = 1 and can be written as a 

trace 

L-\ L-l 

Z LyM = Tr(G LyM A)=Tr( ]T S(M - £ Ti ) J] (nD + (1 - n)E)A). (6) 

{t;=0,1} i=l i=l 

Here 5(x) is the Kronecker S Xt o- In fact, the expression (2) is a conditional probability distribu- 
tion which gives the probability of finding configuration {C} in the stationary state provided 
that a negative particles exists at site L. We will see that all the physical quantities such as 
the speed of particles and the density profile of the positive particles can be written in terms of 
Zl,m- It can easily be checked that a one-dimensional representation of the algebra (3-5) exists 

D = A = E = 1 



2 



for q = 3 [8]. In this case (as can be seen from (2)) all configurations of the system occur 
with equal probabilities P({C}) = z 1 and no correlation exists. It implies that at q = 3 the 
mean-field results are exact. We will see that in this case Zl^m is equal to the total number of 
all possible configurations of the system. Traditionally, one has to find the representations of 
the algebra appears in MPA formalism and use it to calculate both the partition function and 
physical quantities. But in this paper we show that (see Appendix) the expression (6) can be 
calculated directly using the algebra (3-5) and has a closed form 



Zl,m 



i=0 



3)(i)* + i 



C 



M-i 
L-i-2 



(7) 



in which Cf = j\^zjy_ is the binomial coefficient. At q = 3 we have 



M 



ry \ " n M-i n M 

Z L ,M = 2^ C L-i-2 = C L-1 
i=0 

which is nothing except the total number of possible configurations of the model. 



3 Exact Physical Quantities 

In the stationary state, the speed of the positive particles in the reference frame of the lattice 
is found to be (see Appendix) 



V + 



(L - M)Z l _i,m-i + Cf-2 1 
MZ L m 



(8) 
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Figure 1: The speed of the positive particles V+ as a function of q for p = 0.25 (Plot of the 
Eq.(8)). The phase transition takes place at q c = 0.5. 
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Figure 2: The speed of the negative particle V_ as a function of q for p = 0.25 (Plot of the 
Eq.(9)). The phase transition takes place at q c = 0.5. 



Similarly, the speed of the negative particle in the reference frame of the lattice can be obtained 
(see Appendix) 

y = Zl - 1 M-1 + C L-2 + C L-2 ,qn 

Zl,m 

In figures 1 and 2 we have plotted V+ and V-, respectively, computed from the exact expressions 
(8) and (9) as a function of q for L = 500 and M = 125. As can be seen the behaviour of 
these functions changes near q c = 2p. However, we should take the limit (L, M — > oo and 
= p being fixed) to find the exact transition point. We use the steepest descent method 
for computing the transition point. At this limit the asymptotic behaviours of (8) and (9) are 
given by 

v+ =r^ if2 ^ no) 

+ \ I-P if2p<g ' 1 j 

L ^9-2 q-2p 

Thus, the transition point is exactly found to be q c = 2p. Now it is obvious that for q > 2 no 
phase transition can take place because of the restriction on the density p < 1. 
To study the nature of these phases we can compute the density profile of the positive particles 
on the ring. The average density of the positive particles n(x) at the distance x from the 
negative particle is equal to the probability of finding a positive particle at site L — x — 1. Using 
the same procedure proposed in the Appendix, it can be shown that 

, , Zl-im-i ~ g^ggHj r«c£-£$_ 2 }e(M > x) 

n(x) = ^ '- (12) 
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Figure 3: The density profile of the positive particles n(x) for L = 500, M = 200 and q = 0.95 
(Plot of the Eq.(12)). The negative particle is at x = 0. 



where 



G(y > x) 



1 if y > x 
otherwise. 



(13) 



As can be seen from (12) the density profile of positively charged particles is completely flat at 
q = 3. In what follows we study the density profile (12) in the thermodynamic limit. 



3.1 Power-law phase; 2p < q 

In this phase, as can be seen from figure 3, the negative particle has a short-range effect on 
the system. The density profile starts from n(0) just in front of the negative particle and 
decreases exponentially to its bulk value p. The asymptotic behaviour of formula (12) in the 
limit (L, M — ► oo and = p being fixed) is given by 



n(x) = p(l 



(l-p)(q-2)(q-3) 



(q-2)(q-2p)-p{q-2p)-p{q-2Y 



(14) 



in which £ 



In ^| 1 specifies the disturbance due to the existence of the negative particle. 



3.2 Jammed phase; 2p > q 

In this domain, the negative particle provokes a macroscopic shock in the system, as can be seen 
in the figure 4. A high density region puigh = 1 extending from x = to x = xo is separated 
by a sharp interface from a region of low density piow = f extending from x = xq to x = L — 1. 
The position of the shock xq can be obtained from the thermodynamic behaviour of (12) in this 
region 

-2p 



xq = L- 



q-2 
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Figure 4: The density profile of the positive particles n(x) for L = 500, M = 200 and q = 0.45 
(Plot of the Eq.(12)). The negative particle is at x = 0. 



4 Comparison and concluding remarks 

In this paper we investigated a recently introduced model [8,9] in the special limit of infinite 
system size and infinite number of positively charged particles in the presence of one negative 
particle. In this limit, the model shows both phase transition and spatial condensation of par- 
ticles in one-dimension. The phase transition takes place at q c = 2p, between a power-law and 
a jammed phase. The speed of the positive and the negative particles are not differentiable at 
the transition point q c = 2p. In the neutral case, when the density of the positive and negative 
particles are equal, the model shows a first order phase transition at q c = j^fj-, where the 
density of particles as a function of fugacity has a finite jump [8,10]. Very recently it has been 
shown that the transition point can be determined, with great accuracy, using Yang-Lee theory 
in the non-equilibrium context [10]. In [10] the author presents a numerical study, using the 
grand canonical ensemble, but in contrast to the original Yang-Lee paper does not give a proof 
of the existence or nonexistence of the phase transition. Although it was claimed that neither 
in the frame- work of a grand canonical nor canonical ensemble such transition exists [11]. Our 
exact calculations admit that in the limit (L, M — > oo and = p being fixed) the phase 
transition exists and in comparison to the neutral case [8] the transition point is shifted. In 
the neutral case the current of both positive and negative particles is zero for q < 1 where the 
particle segregation exists. Here, as can be seen from Figures 1 and 2, the speed of particles 
never drops to zero for q > 0. In our model the presence of a single negative particle produces 
a macroscopic shock for q < q c . In the power-law phase q > q c the density profile of positive 
particles has an exponential behaviour with correlation length £ which diverges at the transition 
point. This implies a second order phase transition at this point. 
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The next step could be analytical study of this model with unequal but finite densities of positive 
p + and negative particles p-. An approach can be considering the model with open boundaries 
where the injection and extraction rates control the bulk density. Work in this direction is in 
progress [12]. 

Appendix: Formulas for ^l,m 5 V + and V- 

In the following we show that the normalization factor Zl : m can be calculated directly using 
the algebra (3-5). First, we note that one can split (6) in two terms as 

Z l ,m = Tr(G L - 1>M EA) + Tr(G L ^ M -iDA). (15) 

Now using (4) and (5) it is not difficult to see that (15) can be written as 

Z l ,m = Cll 2 Tr(E L ' M - l A) + Tr{G L ^ M -iDA). (16) 

This can be done by noting that the first term in (15) can be written as 

L-2 L-2 

Tr{G L - ljM EA) = Tr( ]T 5(M - ]T n) J] ( n D + (1 - n)E)EA) = 
{T i= 0,l} i=l i=l 

L-2 L-2 

S(M-J2 n)Tr( J] (nD + (1 - n)E)EA) = 

{n=0,l} i=l i=l 

L-2 

6{M-J2 n)Tr(E L ~ M ~ 2 EA) 

{r i= 0,l} »=1 

which gives the first term in (16). One can expand the second term in (16) repeatedly as follows 

Tr{G L - ljM -iDA) = Tr(G L -2,M-iEDA)+Tr(G L - 2 ,M-2DDA) = 

C^Tr(E L ~ M ~ 2 DA) + Tr{G L . 2M -2D 2 A) 
to obtain the following compact formula for the partition function 

M 

Z LM = J2Tr(E L - M - 1 D i A)C^_Z 2 . (17) 

i=0 

The expression (17) can easily be calculated using the algebra (3-5) 

fi = Tr(E L ~ M ~ 1 D i A) = Tr{E L - M - l D i - l \AD + A + D)) = U?fi-i + Tr(E L ~ M ~ 1 )). 



Solving the above difference equation gives 

3)(|) 
q-2 



[<l :i)( " Y f -Tr^-^" 1 ) , i = 0...M. (18) 



7 



Since Zl^m enters to both denominator and numerator of physical quantities, the same common 
factor Tr(E L ~ M ~ r ) cancels from all formulas so the results are independent of its value. For 
simplicity we set Tr(E ) = 1. 

The current of the positive particles is defined as 

J+ = P{n = 1, T l+1 = 0) + qP(n = 1, T l+1 = 2) - P(r t+1 = l,n = 2) (19) 

in which P{ji = m,Tj = n) is the probability of finding the system in configuration {C} 
provided that a particle of kind m is at site i and a particle of kind n at site j. Now we can 
write (19) in term of the conditional probability (2) as follows 

J+ = P{n = 1, r i+ i = 0) + qP{n = l|r m = 2)P(r i+1 = 2) - P(r i+1 = l\n = 2)P(n = 2) = 

p , 1 ^ , 1 f Tr(G L ^ M -iDA) Tr(G l-i,m-\AD) 

P(n = i,n +1 = 0) + I ( q Tt{GlmA) ) 

in which we have used P(n = 2) = \. Using the algebra (3-5) we obtain 

lTr(G L - 2tM -i(A + D)) 



J + = P(n = l,T i+1 = 0) + 



L Tr(G L>M A) 



VI 1 n\ , 1 Z L-1,M~1 + C L-2 /on ^ 
P(T i = l,T i+1 =0) + - (20) 

in which we have Tr(Gi-2,M-iA) = Zl-i,m-i and 

L-2 L-2 

Tr(G L - 2M -iD) = Yl S(M-1-J2 n)Tr( ]J ( n D + (1 - n)E)D) = 
{n=Q,l} »=1 i=l 

E 5(M-l-X:r l )Tr(^-^ 2 ) = Cf_- 2 1 

{Ti=0,l} »=1 

The first term in (20) can be written as 

L 

P(n = 1, t 1+1 = 0) = ]T P(r, = 1, r i+ i = 0|r fc = 2)P(r fc = 2) = 

fc=l,fc^i,i+l 

1 v^ 1 P /_ _ i _ - o\ - 1 V* Tr ( G j,P DEG L-j-i,M-P A ) _ 
L fc=2 L fc=2 P =o J ^C-l,mA) 

(L — M — l)Z L ^ M -i 
LZl,m 

Putting the last term in (20) gives 

(L - M)Z L ^ M -i + Cf^ 1 



J+ = 



LZl,m 
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Now we can calculate the speed of the positive particles 

rr _ L _ (L ~ M)Z L _l jM _l + Cfjt 

V+ - M - MZ^ (21) 

which is exactly the expression (8). The speed of the negative particle is also defined as follows 

F_ = P( n = 0|r m = 2) + qP{r % = l\r i+1 = 2) - P(r i+l = l|r< = 2). 

Now using the expression (2) for the conditional probability and the algebra (3-5) we obtain 

_ Tr{G L - 1M EA) Tr(G L - liM DA) Tr(G L ^ M AD) 
Tr(G LM A) +<1 Tr{G LM A) Tr(G LM A) 

Tr(G L -i }M EA) Zl-^m-x + C^ 1 _ Z L _i >M -i + Cf_ 2 + C^ 1 
Tr(G L>M A) + Z LM ~ Z L . M ■ 1 ' 
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